Direct Measurement of the Spatial-Spectral Structure of Waveguided Parametric 

Down-Conversion 
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We present a study of the propagation of higher-order spatial modes in a waveguided parametric 
down-conversion photon pair source. Observing the multimode photon pair spectrum from a pe- 
riodically poled KTiOP04 waveguide allowed us to isolate individual spatial modes through their 
distinctive spectral properties. We have measured directly the spatial distribution of each mode of 
the photon pairs, confirming the findings of our waveguide model, and demonstrated by coincidence 
measurements that the total parity of the modes is conserved in the nonlinear interaction. Further- 
more, we show that we can combine the advantages of a waveguide source with the potential to 
generate spatially entangled photon pairs as in bulk crystal down-converters. 



The prevalence of parametric down-conversion (PDC) 
as a source of photon pairs is due not only to the high 
quahty of the photons produced but also its experimental 
simplicity relative to other methods of generating single 
photons. Its ubiquity may lead one to believe that PDC 
is a technique with little scope for improvement. How- 
ever, bulk-crystal down-conversion sources suffer from a 
significant drawback: the photon pairs are emitted in 
a cone-shaped pattern making efficient collection diffi- 
cult. This limits both the absolute count rate for a given 
pump power (stimulating the purchase of ever larger and 
more expensive laser systems) and, more importantly, the 
heralding efficiency of any bulk-crystal PDC source. 

By confining photon pair generation to a channel 
waveguide in a nonlinear optical material one can restrict 
the down-converted light to a well-defined set of spatial 
modes rather than allowing the emission to propagate 
at the natural phase-matching angles. This provides a 
straightforward method of controlling the messy spatial 
emission pattern from bulk-crystal down-converters and 
increases the down-conversion collection rate significantly 
[U, 0, H, 0j @|- However, as the wavelength of the down- 
converted pairs is approximately twice that of the pump, 
the light inevitably propagates not only in the funda- 
mental waveguide mode but also in several higher-order 
spatial modes [1, 0, [l] ■ Because of the coupling between 
the spatial and spectral properties of the photon pairs 
imposed by phase-matching, these higher-order waveg- 
uide modes have a significant impact on the spectrum of 
the photon pairs and can markedly degrade source per- 
formance Optimal source design requires that we 
both understand and control the interaction of higher- 
order modes in waveguides [llj. On the other hand, 
the multimode spatial structure in the down-converted 
beams offers new opportunities for advanced quantum 
state preparation ■ Recent experiments utilize entan- 
glement of the orbital angular momentum of the pho- 
ton pairs prepared by bulk crystal PDC, but they rely 
on heavy filter operations by means of holographic state 
selection [ll, 

M,- Waveguided PDC can directly pro- 
vide entangled higher-order spatial modes in analogy to 
orbital angular momentum (0AM) entangled modes of 
photon pairs generated in bulk PDC experiments. 0AM 



modes and their applications have been extensively stud- 
ied recently [H, [l^, [13, [111 with a view to to accessing 
higher-dimensional Hilbert spaces via the generation of 
hypcrcntangled photon pairs [H HO, HH, [12., IS, M ■ 

In this Letter we report the first direct observation 
of photon pairs generated in higher-order spatial modes 
by waveguided parametric down-conversion. We assign 
specific mode labels to each process by applying a nu- 
merical model and confirm parity conservation between 
the interacting mode triplets. Furthermore, we show that 
our source can generate spatially entangled two-photon 
states, while retaining the virtues of a waveguided device. 

In general, down-converted photon pairs from waveg- 
uides are entangled in both frequency and spatial mode. 
Because of the dependence of the mode propagation con- 
stants on wavelength, spatial mode and spectrum are 
linked through the phase-matching conditions. Hence en- 
tanglement in these degrees of freedom cannot usually be 
separated A key property of our source presented 

in this paper is its particular modal dispersion inside the 
waveguide which fulfills all the requirements to generate 
Bell-states in the spatial domain. By spectrally filter- 
ing the down-converted beams this source allows for the 
generation of photon pairs whose spatial entanglement is 
separated from the spectral domain. Hence hyperentan- 
gled photon pairs are emitted: 
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with iV'.s'^'') and \4>f^') denoting the spectral properties 
and I 'f) Bell denoting a Bell state for higher-order spatial 
modes. Equation ([1]) is derived by applying a spectral 
Schmidt decomposition and introducing broadband fre- 



qucncy modes IV's'^^): highlight the decoupling of 

the spectral and spatial degrees of freedom [2^. 

Our source is a 10 mm z-cut pcriodically-polcd 
KTiOP04 (PPKTP) waveguide from AdvR with a nom- 
inal poling period of 8.72 /im, pumped by a spatially fil- 
tered pulsed diode laser at 403.3 nm with a bandwidth 
of 0.8 nm. Input coupling was through a 20x micro- 
scope objective and the pump was observed to be mainly 
(though not exclusively) in the fundamental mode of the 
waveguide. Output coupling was by an aspheric lens with 



a focal length of 6.24 mm, set to image the output face 
of the waveguide to a plane about 800 mm away. The 
type-II phase-matching conditions ensured that we ob- 
tained almost degenerate photon pairs, with the horizon- 
tally polarized pump (?;-polarizcd in the crystal basis) 
yielding signal and idler with horizontal (y) and verti- 
cal (z) polarizations respectively. After the crystal the 
pump was removed with long-pass filters and the signal 
and idler photons were separated at a polarizing beam- 
splitter (PBS). 

Initially, signal and idler beams were coupled into two 
multimode fibers attached to a spectrometer with single- 
photon sensitivity. The multimode fibers allowed us to 
monitor simultaneously a range of spatial modes gener- 
ated in the waveguide. The spectra for signal and idler 
are shown in Fig.[T] The spectral signatures of several 
spatial modes are clearly present: five individual peaks 
can be identified in both spectra. The peaks in the sig- 
nal arm are paired with those in the idler through energy 
conservation. Each pair of peaks corresponds to a par- 
ticular spatial mode set of pump, signal, and idler. Us- 
ing the single-photon spectra, two sets of spectral filters 
(one for the signal with central wavelengths of 808 nm 
and 830 nm and bandwidths of 3 nm and the other for 
the idler with central wavelengths 810, 830, and 860 nm 
and bandwidths of 10 nm) were calibrated such that each 
spectral peak could be individually selected by inserting 
and angle tuning a particular filter. 

In order to assign specific mode labels to peaks A to 
E we developed a model of down-conversion in a step- 
index waveguide with a rectangular profile, bounded by 
a uniform dielectric on three sides and by air at the 
fourth. Although the production method of the waveg- 
uides results in a graded index distribution orthogonal 
to the air interface [J our model yielded a simplified, 
semi-analytic solution which has been proven as suffi- 
cient to describe the experimental results [ll| . Although 
a more precise model would cause slight alterations to the 
predicted peak heights and spatial mode distributions, 
the more salient features — the central wavelengths of 
the spectral peaks — would remain virtually unchanged. 
By adjusting the index contrast, the waveguide dimen- 
sions, and the poling period as free parameters we fitted 
the calculated spectra to the measured marginal spec- 
tra of signal and idler (Fig. [2]). Note that the spectral 
widths of the signal and idler marginal distributions are 
set by the overlap of the pump bandwidth with the modal 
phase-matching functions of the crystal. Here the rel- 
atively broadband pump results in photon pairs with 
wider bandwidths than in similar experiments utilizing 
continuous- wave lasers Q. From the theoretical model- 
ing we identified each of the mode triplets listed in Ta- 
blelH they are labeled with the number of nodes in the 
horizontal and vertical directions (x, y) respectively. 

The principal mode pair A, is the result of the interac- 
tion between the fundamental modes of all three fields. 
This triplet has the most widely separated spectral com- 
ponents and was fitted by adjusting the poling period in 
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FIG. 1: (a) Experimental setup. LD, laser diode; SF, spa- 
tial filter; HWP, half-wave plate; MO, microscope objective; 
KTP, PPKTP waveguide; AL, aspheric lens; LWP, long- 
wave pass filter; PBS, polarizing beamsplitter; Fs, Fi, inter- 
ference filters; FM, flipper mirrors; MMF, multimode fiber; 
APD, avalanche photodiode; GS, grating spectrometer, (b) 
Background-subtracted spectra with correlated peaks labeled 
A to E. Colored lines show filtered spectral peaks (heights 
adjusted for ease of viewing). 



the model. This effective poling period of 8.92 fj,m serves 
as a global correction to allow for the difference between 
the waveguide in the lab and the empirical Sellmeier 
equations p6j . The remaining free parameters were ad- 
justed to reproduce the observed marginal spectra of the 
photon pairs. With waveguide dimensions of 4.1 x 9.3 /im 
and an index contrast of 0.008 we obtained a very good 
agreement between theory and experiment (see Fig. 
These dimensions were verified under an optical micro- 
scope. 

According to our model, mode pairs A, C, and E orig- 
inate from the (0, 0) component of the pump (sec Table 
H]). E stems from photon pairs generated both in modes 
(1,0) and (0,2) with overlapping spectral distributions. 
As a result of this frequency degeneracy and the coher- 
ence of the PDC process, the signal and idler pairs in E 
are entangled in spatial mode. Further mode pairs oc- 
cur at B and D pumped by the fraction of the pump 
intensity in the (0, 1) mode (37.5%). In these cases, the 
signal and idler are in different — though still parity 
conserving — modes. Peaks B and D each consist of two 
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(0,0)p 


-> (0,0). + (0,0), 


B 


(0,l)p 


-> (0,0). + (0,1), and (0,1), + (0,0), 


C 


(0,0)p 


-> (0, 1). + (0, 1), 


D 


(0,l)p 


^ (0,1). + (0,2), and (0,2). + (0,1), 


E 


(0,0)p 


(1,0). + (1,0), and (0,2). + (0,2), 



TABLE I: Processes giving rise to the five observed mode 
spectra. 
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FIG. 2: Comparison of modeled spectra with measured data. 
Labels A to E correspond to labels in Fig.[T] 



down-conversion processes with almost identical spectra 
each entangled in spatial mode. The two pairs of peaks in 
the modeled spectra not mirrored in the measured data 
arc from down-conversion events into higher order spa- 
tial modes [up to (1,2)]. These modes couple poorly into 
the collection fibers and hence are not seen in the data. 
The discrepancies in peak height between theory and ex- 
periment in Fig. [2] stem from our rectangular waveguide 
model and the falling collection efficiencies for higher- 
order spatial modes. 

Next the high-sensitivity CCD camera was removed 
from the spectrometer and placed in the image plane of 
the / = 6.24mm asphcric to measure the spatial intensity 
distribution of each mode. Both output beams from the 
waveguide were directed simultaneously to separate areas 
of the sensor yielding magnified images (approximately 
130 x) of the spatial modes of both signal and idler in 
the waveguide. Fig. [3] shows the characteristic distribu- 
tions of individual spatial modes, recorded by tuning the 
spectral filters to pick out spatial modes through their 
unique spectra. This demonstrates the strong correla- 
tion between the spatial and spectral degrees of freedom 
in this system. 

In these measurements, a high level of background 
was present from the long-lived, unphasematched fluo- 
rescence emitted by the waveguide that could not be 
removed by time-gating due to the slow speed of the 
camera. Instead the fluorescence level was measured for 
each spatial mode by rotating the pump polarization to 



vertical, hence removing any phase-matching. However, 
this background could not be subtracted directly as the 
fluorescence was higher for a vertically polarized pump. 
Therefore, auxiliary measurements for both polarizations 
were made with an unpoled waveguide in which no phase- 
matched processes could take place. The ratio between 
these fluorescence signals allowed us to introduce a cor- 
rection for the scaling of the background in the PDC 
spatial mode images recorded with the poled waveguide. 
Subtracting this adapted background from the spatial 
mode images yielded a realistic measurement of the true 
distribution of the PDC in the various spatial modes as 
shown in Fig.[3l 

It is apparent from Fig. [3] that each of the five mode 
pairs A to E has its own characteristic spatial intensity 
distribution, with peaks B to E exhibiting obvious signs 
of higher-order mode propagation. All of the recorded 
spatial distributions agree very well with those found in 
the spectral degree of freedom through the model as listed 
in TablelH A and C are pure (0, 0) and (0, 1) respectively; 
B is a sum of (0,0) and (0, 1) where two processes over- 
lap spectrally; D is also a sum of two processes, (0, 1) 
and (0,2); E, the superposition of the (1,0) and (0,2) 
modes, is the only case to show a higher-order mode in 
the horizontal direction. Furthermore one can see the 
deviation of the waveguide from a rectangular structure: 
the fundamental mode sits at the top of the guide close to 
the air boundary, while higher-order modes spread down 
into the chip where there is an exponential decay in the 
refractive index contrast not present in our model. 

Finally, we performed a coincidence measurement be- 
tween the different spatial modes. With both beams once 
again coupled into the multimode fibers the photons were 
sent to two silicon avalanche photodiodes (APD). The 
time-gated single count rate of each APD was monitored 
along with the rate of coincidence counts between the two 
as the filters were set to select every combination of the 
five spatial modes in both the signal and idler arms. The 
results are shown in Fig. [3] with the background of acci- 
dental coincidences — calculated as the product of the 
singles rates divided by the laser repetition rate (1 MHz) 
— subtracted from the coincidence rates. The presence of 
only diagonal elements in the corrected coincidence rates 
demonstrates the strict correlation between the spatial 
modes: if the signal photon is emitted into a particular 
spatial mode then the idler will always be found in the 
corresponding mode. This confirms the requirement for 
parity conservation between the three interacting modes. 

Our measurements demonstrate that the generation of 
higher-order spatial mode entanglement can be easily ac- 
complished in waveguided PDC. For our source this can 
be achieved by filtering processes B, D, or E and post- 
selecting on successful coincidence events. For example, 
by choosing only process B we find 



4 



(a) A 




(b) 




Signal mode 



Idler mode 



FIG. 3: (a) Plots of the background-subtracted spatial mode distributions for signal (top) and idler (bottom) arms along with 
their marginal distributions. The air interface is towards the top of each frame, and the decay of the waveguide structure into 
the chip can clearly be observed towards the bottom, (b) Coincidence count rate between different pairs of spatial modes. 
Labels A to E correspond to labels in Fig.[T] 



|^)B = i?'^A,|V'W,#-)) 



(|(0,1)„(0,0),) + |(0,0)„(0,1),)). (2) 



Similarly, filtering peaks D and E yields the Bell states 
1^'"'') and l^"*") respectively. 

In conclusion, we have directly imaged spectrally- 
resolved spatial modes of PDC in a PPKTP waveguide. 
We have identified the individual spatial mode contribu- 
tions and demonstrated that our model accurately repro- 



duces the photon pair spectra. This shows that waveg- 
uided PDC sources have the potential to be used as bright 
sources of photon pairs entangled in spatial mode. These 
photon pairs may have many applications from testing 
the Bell inequalities in the spatial domain to distributing 
Bell pairs over free space links for quantum key distribu- 
tion applications. 
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